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500 million tons of coal fly ash are produced worldwide every year with only 16% of the total amount 
utilized. Therefore, potential applications using fly ash have both environmental and industrial interests. 
Unburned carbon concentration measurements are fundamental to effective fly ash applications. Current 
on-line measurement accuracies are strongly affected by the mineral content and coal rank. This paper 
describes a char/ash particle cluster spectral emittance method for unburned carbon concentration 
measurements. The char/ash particle cluster spectral emittance is predicted theoretically here for various 
unburned carbon concentrations to show that the measurements are sensitive to unburned carbon 
concentration but insensitive to the mineral content and coal rank at short wavelengths. The results show 
that the char/ash particle cluster spectral emittance method is a novel and promising route for unburned 
carbon concentration on-line measurements without being influenced by mineral content or coal rank 
effects. 



The world currently produces 500 million tons of coal fly ash each year 1 . This fly ash is quite dangerous 
because of its high toxic element content and which can affect people's health and cause genetic mutations 2 " 7 . 
Therefore, fly ash is collected and either utilized or disposed of in landfills at significant costs 8 ' 9 . Currently 
only 16% of the total fly ash is utilized worldwide 1 which means that significant disposal costs are charged to the 
power companies and, thus, to consumers. Therefore, research on potential applications to use fly ash has both 
environmental and industrial interests. 

The high specific surface area makes fly ash a low cost adsorbent for coal fired flue gas and waste water 
treatment 1 . Intrinsic fly ash and activated carbon and zeolite generated from fly ash have been used to capture 
C0 2 10 " 17 and to adsorb SO x 18 , NO x 18 and mercury 19 " 25 in flue gas and to remove toxic metals and organic 
compounds 26 " 29 from waste water. The geotechnical properties (specific gravity, permeability and consolidation), 
make fly ash a promising construction material 1 . Investigations have shown that the unburned carbon compon- 
ent in fly ash is important for its applications 110 ' 18,21 ' 22,25 . Fly ash must have a low unburned carbon content when 
used as a construction material in the cement industry since the unburned carbon undesirably adsorbs the air 
entraining agents added to prevent cement crack formation 10 . The unburned carbon concentrations in fly ash 
have increased dramatically since utility boilers have been using low NO x burners to meet the Clean Air Act NO x 
emission requirements. As a result the unburned carbon concentrations in fly ash range from 2% to 20%, which 
restricts fly ash use as a construction material. However, the higher unburned carbon concentrations make fly ash 
a more efficient adsorbent for flue gas and waste water treatment 18 ' 21 ' 22,25 since the unburned carbon contributes 
much of the surface area to fly ash. Therefore, unburned carbon concentration measurements are of fundamental 
importance in fly ash applications. 

The standard unburned carbon concentration measurement method is the loss-on-ignition test. Sample fly ash 
is removed from the exhaust, and then oxidized at 850°C for several hours in a muffle furnace. The weight loss is 
the unburned carbon weight. However, the analytical process takes up to 24 hours so it cannot be used on-line as a 
real time monitor of the unburned carbon concentration and the boiler working status. The photo-acoustic effect 
has also been used to measure unburned carbon concentrations 30,31 . However, the fly ash sample has to be ground 
prior to measurements to reduce the carbon particle size effect 30 , so this method cannot be used on-line. On-line 
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methods include microwave absorption 32,33 , optical dispersion- 
reflection 34,35 , laser-induced breakdown spectroscopy 36,37 and infra- 
red emissions 38 . The measurement accuracies of these methods are 
strongly affected by the mineral content and the coal rank 36,37,39 . The 
relative errors can be as high as 30% for the microwave absorption 
method 32 , 50% for the optical dispersion -reflection method 34 , and 
40% for the laser-induced breakdown spectroscopy method 37 . 

This paper describes a char/ash particle cluster spectral emittance 
method to measure the unburned carbon concentration. The char/ 
ash particle cluster spectral emittance is shown theoretically to be 
sensitive to the unburned carbon concentrations and insensitive to 
the mineral content and coal rank at short wavelengths. Thus, the 
char/ash particle cluster spectral emittance method is a novel and 
promising route for on-line unburned carbon concentration mea- 
surements that are insensitive to the mineral content and coal rank 
effects. 

Results 

Unburned carbon concentration effect. Figure 1(a) shows the 
spectral emittance of char/ash particle clusters of a typical Chinese 
bituminous coal, Shenhua coal, produced in Shanxi Province, with 
unburned carbon concentrations, C UB , ranging from 0 (pure fly ash 
cluster) to 20%. The char and ash particle sauter diameters were both 
9 um, the particle volume fraction was 0.2 and the mineral content 
are summarized in Table l 40 . 

The spectral emittance is low at wavelengths below 4 um, 
increases sharply at wavelengths between 4 and 8 um, and 
approaches unity at 8 urn for the pure fly ash cluster and the char/ 
ash clusters with C UB below 5%. The spectral emittance is almost 
constant within the whole wavelength range for the char/ash particle 
clusters with C UB above 10%. This evident grey nature agrees with 
spectral emittance experimental data for devolatilized carbon-rich 
char particle clusters 41 . The spectral emittance is close to unity and 
does not change with C UB at wavelengths above 7 um. However, the 
spectral emittance is strongly sensitive to C UB at wavelengths below 
4 um, which implies that C UB can be determined by measuring the 
char/ash particle cluster spectral emittance at a low wavelength. The 
char/ash particle cluster spectral emittance at 4 um is plotted versus 
C UB in figure 1(b). The slope is higher for C UB below 10%; thus, the 
char/ash particle cluster spectral emittance method tends to be more 



Table 1 | Mineral content of Shenhua and Yankuang coal fly ash 

Si0 2 /% Al 2 0 3 /% CaO/% Fe 2 0 3 /% 

Shenhua 56.83 26.17 10.20 6.98 

Yankuang 44.50 35.00 13.02 7.48 



accurate for measuring the unburned carbon concentration in fly ash 
with lower C UB . 

Mineral content and coal rank effects. Current on-line measure- 
ment methods are strongly affected by the mineral content and coal 
rank 32-38 . These effects are investigated here for the char/ash particle 
cluster spectral emittance method. Figure 2(a) compares the spectral 
emittance at 4 um for char/ash particle clusters of Shenhua coal and 
another typical Chinese bituminous coal, Yankuang coal, produced 
in Shandong Province, which have different mineral content as 
summarized in Table l 42 . The average relative difference of the 
spectral emittance for these two coal char/ash clusters is 0.4%, 
which shows that the mineral content has little effect on the char/ 
ash particle cluster spectral emittance and, thus, on the unburned 
carbon concentration measurements using the char/ash particle 
cluster spectral emittance method. 

For low rank coals such as sub-bituminous coal and lignite, the fly 
ash has more calcium oxide and less silica than high rank coals 1 . The 
unburned carbon content also has larger particles and higher vitri- 
nite maceral components which are more porous than inertinite 
components; thus, lower rank coals have higher porosities 43 . As 
explained above, the mineral content has little effect on the measure- 
ments. Since the carbon-rich extracts separated from the fly ash 
mineral content for low rank coals also typically have carbon con- 
centrations of about 75 % 44 ' 45 , the carbon concentration of a single 
char particle and the number of char and ash particles remained 
unchanged. The porosity of a single char particle, however, increased 
from 50% to 75%, and the particle diameter of a single char particle 
increased from 9 um to 11 um accordingly for the Shenhua coal 
char/ash clusters to investigate the effect of coal rank. The average 
relative difference of the spectral emittance before and after changing 
the porosity and particle diameter is 2.5% as illustrated in figure 2(b), 
which shows that the coal rank has little effect on the char/ash par- 
ticle cluster spectral emittance and, thus, on the unburned carbon 
concentration measurements. 

Implications. The char/ash particle cluster spectral emittance 
depends on the optical depth, t, the phase function, 0(in.[i'), and 
the scattering albedo, w, as shown in equation (1). Char/ash particle 
clusters are opaque with a thickness on the order of 10 2 um 3 , while 
the thickness is generally on the order of 10 4 um for spectral 
emittance measurements, which means that t has little effect on 
the char/ash particle cluster spectral emittance. However, the 
spectral emittance was affected by t for spectral emittance 
measurements of suspended semi-transparent char/ash particle 
clouds 38 . The mineral content and coal rank then affect t and thus, 
the spectral emittance and unburned carbon concentration 
measurements. Therefore, the optically opaque nature of char/ash 
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Figure 2 | Effects of (a) mineral content and (b) coal rank (the two porosities represent different coal ranks) on the spectral emittance. 



particle clusters is the first reason that the mineral content and coal 
rank have little effect on the unburned carbon concentration 
measurements. Char and ash are both strongly forward scattering 
particles 46 regardless of the mineral content or coal ranks. Therefore, 
the mineral content and coal rank have little effect on &{fi> fi') and, 
thus, on the spectral emittance and unburned carbon concentration 
measurements. 

Since the optical depth and phase function have little effect, the 
char/ash particle cluster spectral emittance is only affected by the 
scattering albedo, which is further seen in the similar trends of the £h 
~ k and (1 — co) ~ k curves plotted in figure 3(a). 1 — co charac- 
terizes the ratio of the absorption efficiency to the extinction effi- 
ciency. Therefore, measurements of the char/ash particle cluster 
spectral emittance are actually relative measurements, unlike the 
absolute measurements of the transmission power in the microwave 
absorption method 32 ' 33 , the reflected light intensity in the optical 
dispersion-reflection method 34,35 , and the atomic emission signal in 
the laser-induced breakdown spectroscopy method 36 ' 37 . The similar 
trends in the k ~ k and (1 — co) ~ k curves plotted in figure 3(b) 
show that the absorption index, k, is the dominant effect to deter- 
mine the ratio of the absorption efficiency to the extinction effi- 
ciency, and, thus, the cluster spectral emittance. Given that the 
absorption index of carbon is three orders of magnitude higher than 
that of pure ash at wavelengths below 4 um, the char/ash particle 
cluster spectral emittance at short wavelengths is sensitive to the 
unburned carbon concentrations and insensitive to the mineral con- 
tent and coal rank. The fly ash cluster spectral emittance has been 
measured on-line by Moore et al. 47 . Therefore, this paper provides 
implications that the char/ash particle cluster spectral emittance 
method is a promising route for unburned carbon concentration 



on-line measurements that is independent of the mineral content 
and coal rank effects. 

Discussion 

This paper describes the char/ash particle cluster spectral emittance 
method for unburned carbon concentration measurements. The 
effects of mineral content and coal rank were investigated to deter- 
mine their influence on the spectral emittance measurements. The 
spectral emittance was found to be sensitive to the unburned carbon 
concentrations and insensitive to mineral content and coal rank at 
short wavelengths due to the optically thick and strong forward 
scattering nature of char/ash particle clusters, and the fact that the 
absorption index of carbon is three orders of magnitude higher than 
that of pure ash at short wavelengths. This paper provides implica- 
tions that the char/ash particle cluster spectral emittance method 
described here is a novel and promising route for on-line unburned 
carbon concentration measurements that is not influenced by min- 
eral content and coal rank effects. 

Methods 

A char/ash particle cluster is modeled as an isothermal, one-dimensional, equivalent 
homogeneous semitransparent absorbing and scattering layer. The radiative transfer 
equation is 48 : 



dx 



:(l- W )/ b -/(T,/i)+|j ®{W')I{T,V'W (1) 



where t is the optical depth, J(i,/i) is the radiation intensity at position x and direction 
H (=cos0, where 6 is the angle between the direction and the positive x direction), I b is 
the radiation intensity emitted by the blackbody at the cluster temperature, 
co ( = Q S /Qe> where Q s is the scattering efficiency and Q e is the extinction efficiency) is 
the scattering albedo and //) is the phase function. 





Figure 3 | (a) Spectral emittance and (b) Absorption index change with scattering albedo. 



SCIENTIFIC REPORTS | 4:4567 | DOI: 1 0.1 038/srep04567 



3 



The radiative transfer equation was solved using the S-12 discrete ordinate 
method 49 with equation (1) then converted to 24 differential equations, one for each 
transfer direction, expressed as: 

p, j / \ 24 

ft-^ = (l-o))I b -I(T,^)+ WJ^^W^) (2) 

The boundary conditions are: 

/(0,/i>0) = 0 (3) 
/( ^ w,m) = l( i w, - (4) 
Combining of equations (2)-(4) gives the hemispherical emittance, e h > as: 

24 

£ w(w,/i>o) 



E Wb 

Fiveland 49 gave details on how to solve equation (2). 

The radiative transfer equation model requires knowledge of the scattering albedo 
and the phase function to calculate the char/ash particle cluster emittance using 
equation (5). Three common cluster radiative property computation methods are the 
multiple sphere T-matrix, discrete dipole approximation (DDA) and generalized 
multiparticle Mie- solution (GMM). The T-matrix method relies on more accurate 
equations for every particle's optical cross- sections than the DDA method, and is 
capable of multi-processor parallelization that is not possible with the GMM 
method 50 . Therefore, the multiple sphere T-matrix method was implemented in the 
open source code by Mackowski and Mishchenko 50 to accurately calculate the cluster 
radiative properties. 

The T-matrix method uses the size parameter x ( = 7iDM, where D is the particle 
diameter and X is the wavelength), optical constants m ( = n + Ki> where n is the 
refractive index and k is the absorption index) and position profiles of each particle in 
the cluster. This algorithm has several steps 50 : (1) Determine each particle's individual 
scattering parameters using Mie theory as expansions of the vector spherical wave 
functions (VSWF). (2) Use the extended boundary condition method (EBCM) to 
calculate the relationships between the expansion coefficients for different particles as 
a set of linear equations. (3) Iteratively solve the resulting equations until the residual 
satisfies the required error limit. (4) Use the scattering from all the particles to 
calculate the radiative parameters of each particle. 

Representative aggregates were generated to provide realistic particle geometric 
distributions in the char/ash clusters for the T-matrix calculations. These aggregates 
are representative of actual char/ash particle clusters in terms of the particle size, 
porosity and unburned carbon concentration. Representative aggregates were first 
generated for high rank bituminous and anthracite coals. The coal rank effect was 
then analyzed in the results and discussion section. 

The fly ash carbon has been found to be distributed in discrete carbon-rich char 
particles 51 , so the representative aggregates included carbon-free ash particles and 
carbon-rich char particles. The ash and char particles have been found to have similar 
size distributions for high rank coals 52 . The effect of the particle size distribution has 




Figure 4 | Spectral emittance of various representative aggregates. The 

left inset is the 524 particle aggregate generated using the Bentz method, the 
middle one is a 524 particle aggregate generated using the DLA method, 
and the right one is an 1114 particle aggregate generated using the Bentz 
method. 



been determined by comparing the spectral emittance calculated using Mie theory for 
ash clusters with uniform diameters and with lognormal size distributions 3 . This 
effect was found to cause average relative differences of less than 3%. Therefore, the 
ash and char particles were assumed to have one uniform diameter. 

Hurt et al. 44 ' 45 analyzed carbon-rich extracts from the fly ash mineral content of 
eight types of coal to show that these extracts had a typical carbon concentration of 
75%. In addition, the unburned carbon content has a dominant maceral component 
of inertinite 53 with 50% porosity 43 . Therefore, each char particle was assumed to be a 
50% porosity sphere consisting of 25% ash and 75% carbon. The Maxwell- Garnett 
effective medium theory 3 ' 44 was used to determine the ash and char particle optical 
constants from the chemical compositions. The representative aggregates were gen- 
erated using the open- source code by Bentz 55 . The particles were first randomly 
placed in a cube and then each particle walked randomly until it contacted and 
adhered to another particle. The numbers of ash and char particles were determined 
by the unburned carbon concentrations and particle volume fractions. 

Figure 4 shows the calculated spectral emittance of representative 524 and 1114 
particle aggregates generated using the Bentz method and a 524 particle aggregate 
generated using the diffusion-limited aggregation (DLA) method 56 for a char/ash 
particle cluster with 9 um diameter, 0.2 particle volume fraction and 2% unburned 
carbon concentration. The optical constants were taken from Goodwin and 
Mitchner 57 for the mineral content, Si0 2 , Fe 2 0 3 , A1 2 0 3 , CaO and MgO, and from 
Foster and Howarth 58 for the carbon. These results agree well with each other which 
means that the 524 particle aggregate is representative of the char/ash particle cluster 
and that the different random aggregation algorithms have little effect on the cluster 
radiative properties. 

Further validation of the theoretical method and results is provided in the 
Supplementary Information. 
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